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Abstract 

We determine the volume and mass dependence of scalar and pseudoscalar two-point functions in Nj- 
flavour QCD, in the presence of an isospin chemical potential and at fixed gauge- field topology. We obtain 
these results at second order in the e-expansion of Chiral Perturbation Theory and evaluate all relevant zero- 
mode group integrals analytically. The virtue of working with a non- vanishing chemical potential is that it 
provides the correlation functions with a dependence on both the chiral condensate, E, and the pion decay 
constant, F, already at leading order. Our results may therefore be useful for improving the determination 
of these constants from lattice QCD calculations. As a side product, we rectify an earlier calculation of the 
0{e^) finite-volume correction to the decay constant appearing in the partition function. We also compute 
a generalised partition function which is useful for evaluating U{Nf) group integrals. 
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1 Introduction 



With the possibility of performing unquenched lattice gauge theory calculations at increasingly small quark 
masses, new opportunities for understanding the low-energy properties of QCD are emerging. The theoretical 
framework for such analyses is provided by Chiral Perturbation Theory (ChPT), which describes the long wave- 
length excitations of QCD above the vacuum, in terms of physical parameters known as low-energy constants 
(LECs). While there exist a number of ways to approach the chiral limit, we will consider here the e-regime 
of [1]. This regime is well suited to numerical lattice calculations, since it is reached by lowering the quark 
masses in such a way that the Compton wavelength of the corresponding pseudo-Goldstone bosons is kept much 
larger than the typical linear extent of the finite (simulated) four-volume V. Low energy constants can be 
determined in this regime, by comparing lattice calculations of a given observable with analytical expressions 
for this observable obtained in ChPT. In the present work we will focus on ChPT calculations which allow the 
determination of the two leading LECs, the chiral condensate S and the pion decay constant F, and we refer 
to [2] for a recent review of lattice results for these as well as higher order LECs. 

Two-point correlation functions, of either scalar and pseudoscalar densities, or of vector and axial vector 
currents, are relatively straightfoward quantities to compute with lattice QCD methods. They are also calculable 
in the e-regime of ChPT, and have been studied up to and including second order corrections, in the quenched, 
partially quenched or full theory, with or without fixed gauge-field topology [H HI [5l O [7] . Comparisons of these 
analytical predictions with lattice results have been reported in [8]. At second order, these correlation functions 
couple to both LECs, S and F. However, in scalar and pseudoscalar two-point functions, F only appears at 
sub-leading 0{e^) ~ 1/VV. In vector and axial correlation functions, it is present at leading non-vanishing 
order, but these functions are C'(e^)-suppressed compared to the scalar and pseudoscalar correlators. All of these 
suppression factors are clearly a handicap for determining F from lattice calculations of the above two-point 
functions. 

For that reason, we propose to consider two-point scalar and pseudoscalar correlation functions in the pres- 
ence of a non-zero imaginary isospin chemical potential /i, leaving vector and axial-vector correlation functions 
for future work. The interest of adding a chemical potential is that it couples the decay constant, F, to the zero 
modes in the leading order Lagrangian, through the combination ^^VF^ [S]. This implies that the scalar and 
pseudoscalar correlation functions now depend on F, already at leading order. In turn, the increased sensitivity 
to F should facilitate the extraction of this LEC from lattice calculations. Since all of our results depend only 
on /i^, they hold true for both real and imaginary isospin chemical potentials for any number of flavours Nf. 
For real /i, only the even flavour case has a real action. 

The proposal to use a chemical potential to facilitate the determination of F was originally made in the 
context of spectral correlators [lOl [TT] . In [10] , the spectral two-point density was singled out as particularly 
sensitive to F. More recently, all spectral correlators [12] and individual eigenvalue correlations [13| with 
imaginary isospin chemical potential were computed in chiral Random Two-Matrix Theory, which is equivalent 
to leading order ChPT in the epsilon regime [14j . This setup includes the partially quenched theory, where 
/i-dependent valence quarks propagate in a background of /i-independent sea quarks, allowing for the use of 
standard, zero-density lattice configurations. In fact, a lattice calculation of F has already been performed in 
this framework |15j . 

The remainder of the article is organised as follows. In Sec. [U we derive the C(e^) improved partition 
function and two-point meson correlation functions in terms of expectation values over the pion zero-mode. 
In particular, we correct an earlier calculation ^16] of the one-loop flnite-volume corrections to the partition 
function. In Sec. [3] we present our results for both flavoured and unflavoured scalar and pseudoscalar two-point 
correlation fimctions in TVj-flavour QCD, given in terms of the finite-volume propagator, the mass dependent 
condensate and other generating functions. These general expressions are spelled out most explicitly in the case 
of two flavours. Our conclusions are given in Sec. [H Two appendices. [K\ and [Bl provide technical details for 
group integral identities, including a new result for a generalised partition function. 
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2 The epsilon expansion in ChPT with a chemical potential 



Since imaginary chemical potentials couple to fields like the time component of a constant vector current, they 
can be accounted for in ChPT through the covariant derivative on the pion fields, U{x) G SU{Nf): 



VpU{x) = dpU{x)~iSp^o[C,Uix)] , Vpt/(x)t ^ dpUixY -iSp^o[C,U{xy] 



(2.1) 



where C is a matrix proportional to the chemical potential fi, whose form depends on the charge ^ couples to. 
Moreover, it was shown in [1], for calculations up to and including second order in the e-expansion, that the 
usual leading order chiral Lagrangian is sufficient. Thus, we consider here the Lagrangian: 



£2 = ^Tr [VpUixyVpUix)] - |Tr [M^Uix) + Uix^M] , 



(2.2) 



where M is the mass matrix. We will take M. = diag(mi, . . . ,771^?^) and throughout this section the external 
vector current, C, will be kept general. In Sec. [3] we will specify this current to be the isospin charge matrix, 
C = diag(/iljvi, — /^ljV2) = M^- The simplest and interesting case of Nf = A^i + iV2 = 1 + 1 flavours of equal 
mass mi^2 = w will be displayed in great detail in Sec. [31 As our final results will depend only on fi^, we 
can also rotate back to real chemical potential by fi^ —fi^ at the end of our calculations, for any number of 
flavours Nf = Ni + N2. 

We use the standard parametrisation of the Goldstone flelds for the e-cxpansion: 



U{x) = Uq exp 



(2.3) 



These fields live on the Goldstone group manifold, U{x) £ SU{Nf), whereas the propagating Hermitean fields 
£,{x) = ^(x)^ are members of the corresponding Lie algebra, e.g. for two flavours, Nf = 2, they are given by 
£,{x) = ■^(Ja£,a{x), in terms of the Pauli matrices. In particular, we have split off the zero momentum mode, 
Uq € SU{Nf), explicitly, and will treat its effects exactly. The power counting in the epsilon expansion in the 
presence of a chemical potential is given by [1] : 



V 



M 



do 



ax) 



(2.4) 



All other quantities, such as E and F, are of 0(e"). 

With this power counting, it is straightforward to expand the action, S2 = J d'^x C2 to O(e^). With the 



notation ^2 
• 0{l): 



5(0) + 5(1) + 5(2) _^ 0(g3)^ obtain: 



5(0) = S. 



(0) 
Ua 



c(0) 



4?(S,F) 



[C,ul][C,Uo] 



c(0) 

Oq2 



F^ 
-V—Tt 
4 

i / d'^x Ti-[dp^{x)dpax)] 



V-Tt[M^Uo 



(2.5) 
(2.6) 



Oie): 



S 



(2) 
92 



5(2) = 

1 



s. 



(2) 



12f2 
1 



5(1) = 



S. 



(2) 



-- / d-'xTY 



V2F 



d'^x Tr 
d'^x Tr 



dpax),ax)][dpax),ax)] 



J d^xTT[M^Uoax)'+ax?U^,M] 



d'^x Tr 



[/tcc/oK(x),[c,e(x) 



(2.7) 



(2.8) 



(2.9) 



where the superscripts specify the order in the e-expansion. Note that all terms containing C vanish for C oc IjV/ , 
i.e. pions do not couple to baryon chemical potential. 
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2.1 Partition function at C(e^) and effective couplings 

The partition function for the A^/-flavour chiral theory is given by 



SU{Nf) 

r 

SU{Nf) 



[dHU{x)] exp[— 5] 
dHUocxp[-S^^^i^,F)]Z^iUo) , 



(2.10) 



where 



Z^iUo) ^ J mx)] m exp[4°^(S, F) - S] 



(2.11) 



and where S is the chiral action at an, as of yet, unspecified order. In Eq. (j2.10|) we integrate over the Goldstone 
manifold with Haar measure cIhU and in Eq. (|2ll|) . ,7(0 = {1 - ^piy J d'^xTil^ix)^] + 0(e^)} is the Jacobian 
corresponding to the change of variables of Eq. (|2.3p [1] . 
In this theory, expectation values are given by 



1 

Z 



SU{Nf) 



[dHU{x)] O exp[-5] 



(2.12) 



As we are interested in evaluating observables up to 0{e'^), we need the partition function at that order. Since 
iS^°'(S, F) is 0(1), we must compute Z^{Uo) to 0{e^). Expanding J(^) exp[5^"' (S, F) — S] to second order and 
performing the resulting Gaussian integrals with the propagator. 
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we obtain 



{N] - 1) 



2F^ Nf 
VNf 



Nf 



A(0)Tr[7Wtc/o + C/t7W)] 



Sijdki)A{x - y) 



(2.13) 



A(0)- / d*x[5oA(a;)]^ Tr [C, C/T][C, C/q 



(2.14) 



where M is an overall normalisation factor which does not contribute to the expectation values defined above. 
Now, re-exponentiating the corrections in Eq. (|2.14p . we obtain for the partition function: 



SU{Nf) 



dnUo exph4" (Ee//,i^e//)] +0(e') 



(2.15) 



with E and F replaced in the argument of slf^ by the 1-loop corrected couplings: 



"eff 



F 1- 



F2 



A(0) , 



A(0) - / d^x[doAix)]- 



(2.16) 
(2.17) 



While the correction to the condensate E Eq. (|2.16p has been known a long time [IJ, the correction to F 
was computed only very recently |16j . apart from the second term in Eq. (j2.17p which is new and seems to 
have been omitted in |16j . This term arises from the contribution proportional to {S^^'^)'^ in the computation 
of Z5([/o) to ©(£2). 

In dimensional regularisation, the propagator A(0) is finite and is given by |17j : 



A(0) = -Pi/VV 



(2.18) 
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Moreover [3], 



2VV 



f3i 



J.2 



V 



(2.19) 



where T is the time extent of the box in which the system is enclosed. In Eqs. (|2.18p and (j2.19p . /3i and fcoo 
are numerical constants which depend on the geometry of the box considered [I71[3]. Together with Eqs. (j2.16p 
and (|2.17p . these equations imply: 



'eff 



= E l 



F 1 



NfF^VV , 



Pi 



V 



Coo 



Nf 



2F2/K 



In the particular case of hypercube, i.e. a box with sides T—Li~L2^L'i~V^/^ , fcoo = /3i/2- 



2.2 Two-point correlation functions 

We consider here two-point correlators of the scalar and pseudoscalar quark bilinears, 

So{x) = iij{x)lNfi>ix) , Sbix) = i>{x)tblNfipix) , 
Po{x) = ^l){x)i-1^lNi^{x) , Pb{x) EE ■tp{x)tbij5'i-Nfip{x) , 



(2.20) 



where the ta denote the SU{Nf) generators, and we have normalised Tr[(tQ)2] = i. For Nj = 2 we have 
ta = \<7a in Eq. ()2.20p . the Pauh matrices for a = 1, 2, 3, and \ times the identity for ta=o- 

In the effective theory, the scalar and pseudoscalar densities are most easily obtained by introducing Her- 
mitcan sources, s{x) = Sa(x)ta and p{x) = Pa{x)ta, which have the same spurionic transformation properties as 
in QCD; i.e. through the replacement [T5| : 



M ^ M + s{x) +ip{x) . 



(2.21) 



The e-expansion counting for the sources is thus the same as that of the quark masses, i.e. s{x), p{x) ^ e^. 

Once this replacement is made, the partition function depends on the sources, and the two-point functions 
are obtained, as usual, by taking adequate functional derivatives: 



{Sa{x)Sb{0))uix) 



1 



Z Ssa{x)Ssb{0) 



Z[s,p] 



s—p—0 



(2.22) 



and likewise for the pseudoscalar correlators. Expanding the observables, the action and the Jacobian to ©(e^), 
a calculation analogous to the one performed for the partition function in the preceeding section yields: 



{So{x)Somui. 



^^^^/(Tr[[/o + C/t])^) 



y2 



2F^ 



1 



(Tr {Uo - )u„ - ((Tr[[/o - U^o])\ 

I J l\f \ I Uo 



4 

+ O(e'), 

where all that remains are expectation values with respect to the zero mode, given by 



{O)uo ^ 



1 



dnUoO e "^^o 



Aix) 
(2.23) 

(2.24) 



In the second term of Eq. (|2.23p . since A{x) is ©(e^), the expectation values with respect to Uq can be calculated 

with (E, F) instead of 5^,^ (Ee// ^Feff)- _ 

Along the same lines, one can derive a similar expression for the pseudoscalar correlator: 



(Po(i)fl,(0))i,,,, 



^eff 



(Tr (L/o + t/J 



t^2 



).o-^((Tr[t/o + C/t])^)^^ 



Aix) 
(2.25) 
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Both results Eqs. ([2?23l) and ([2?25l) have the same form as the corresponding expressions at zero chemical 
potential, A* = [I]- However, when inserting the group integral averages over Uq below, they will differ 
explicitly by /i-dependent terms. 

Flavoured two-point functions can be computed exactly in the same way. But there is an important difference 
to the calculation in [J: while for fi ~ and equal quark masses one has {Sa{x)Sb{0))ij(^x^ ~ Sat, this is explicitly 
broken by the chemical potential ^ 7^ (as well as by non-degenerate masses of course). 

For that reason we will only compute the following diagonal sum over flavoured combinations, where we can 
use the SU{Nf) completeness relation J2ai^a)ijita)ki 



A{x) 



(Tr \iU, + Ul)%, ((Tr[[/o + U^])') 

L J iV/ \ / Uo 



4i?2 



+27V| + ^(Tr[(C/o-[/t)2])c/„ 



\{Nj + 2){{Tr[Uo-U^o]f) 



Uo 



(2.26) 



For clarity we have made the sum over flavour indices explicit (we always use summation conventions unless 
otherwise stated). In complete analogy, one obtains for the flavoured pseudoscalar correlation function. 



Uix) 



(Tr \{U, - Ulf])u, - ((Tr[L/o - Ul]f) 

L J iV f \ /Co 



-A(x) 



y2 



^^((Trlt/o-C/t])^^^^ 



^^[N] + 2)[{T,[U, + Ul]f) 



Uo 



-2Nj 



t^2 



— iTr[iUo + U^o) 



lUo 



(2.27) 



Both results agree again with [4] at = 0, when normalised by [N'j — 1). The resulting group averages differ 
however, both by explicitly /x-dependent factors and by a functional change of the mass dependent condensate. 
This is the subject of the next section. 

Before turning to the evaluation of the relevant group integrals, it is useful to make a comment on the 
spacetime dependence of the above correlation functions. For comparisons with lattice QCD calculations, it 
is useful to consider the corresponding zero-momentum correlation functions, which arc functions only of the 
Euclidean time t, i.e. 



T 



Cs{t) = - I d'x{Sa{x)SamU(.) 



and 



T 



Cp{t) = - I d'x{Paix)Pamu(.) , 



(2.28) 



where either a = 0, for the singlet case, or there is an implicit sum over the SU{Nf) adjoint index a, for the 
flavoured case. As above, in Eq. (|2.28p T is the time extent of the box in which the system is enclosed and V/T 
is its spatial volume. For the massless propagtor [T7] . 



1 



/ii(r)^- ld'xA{x) = - 



^ 2' 12 



(2.29) 



with < T < 1 and where t = t/T. Thus, the zero-momentum correlation functions Cs.p(t) are simply obtained 
from the expressions of Eqs. ([123, (P?^ and f^n7\i above, and those of Eqs. ^J^, and 

(j3.25p below, by making the replacement: 



A{x)-^^h,{T) 



(2.30) 
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3 Results for scalar and pseudoscalar correlators at fixed topology 

In order to make the group integrals, which appear in the results of the previous section, tractable analytically, 
we choose to work in sectors of fixed gauge-field topology. This is done by introducing the thcta vacuum angle, 
9, of QCD into the effective theory, through the replacement [T5] 



M + s{x) + ip{x) 



(3.1) 



and Fourier transforming the results with respect to this angle. Thus for instance, the partition function in the 
sector of topology v is given by: 



It is easy to convince oneself that such a transformation has for effect to replace all SU{Nf) group integrals by 
integrals over the much simpler group manifold of U{Nf). Thus, to O(e^): 



dHUo det[[/Q ] exp 



U{Nf 



lvF^ff^,'Tr[^U^rUo] + iFSe//mTr[C/o + U^] 



0)k ^ 



u JU(Nf) 



dHUaOAei[U^]eyiY> 



+ 0{^), (3.3) 



(3.4) 



where we have absorbed the constant M of Eq. (|2.15p into the measure dnUo, as well as the constant term 
\VF^ff^,'TT[T'] from5^°)(Se//,Fe//) 



In the following, we only consider the theory with equal masses, 



(3.5) 



This includes the interesting case of Nf ~ 2 with degenerate up and down quarks. Furthermore, we specify the 
external vector current to be 



C 



/iF , with r = diag(lArj,— 1 



N2 



(3.6) 



such that F^ = Ij^^ and we will also consider more general partition functions, containing higher powers 

Tr[(FC/QFC/o)'°] as generating functionals, for computing all group integrals. 

Correlation functions at fixed topology are obtained by replacing the averages ( O )f7o ^ ( O in Eqs. (|2.23p 

and (HHH). The building blocks that we have to compute in Eqs. ((T^ and ((T^ are ((Tr[;7o])^ + (Tr[f7j;])2\'^ 



(7o' 



(Tr[C/o]Tr[f/J;])^^, and (Tr[;7o2] + Tr[C/o^^])^^. With the help of some U{Nf) group integral identities derived in 
Appendix El these can be expressed in terms of the following known quantities: the partition function and its 
derivatives with respect to mass, chemical potential and an additional external field. The equations that we 
need are: 



(Tr[C/o - U^o])^ 



= -4 



dr]eff S 

Nf S.(7?e//,«e//) , 4i^'^/ 



2 Ar2 



S2 



' Uo Veff S 

for the squared traces, and for single traces of the squared matrices 

t2 



^eff 



(3.7) 

(3.8) 



(Tr {Uo±Wo) 



t^2 



= -4 



Veff 

4a? 



^eff 



+ {2±2)Nf 



4a ff / 1 
+^[Nfy.{Veff,aeff) - (Tr[F])2 + -alj^iMVeff^a^ff) - Nf) 

Veff ^ ^ 



(3.9) 



The generating functionals on the right hand sides are as follows. The mass dependent condensate is given by 

T.^{rjeff,aeff) _ 1 d 



NfOveff'-'^^^"^ = 4(Tr[f/o + C/ot])^.„. 



(3.10) 
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It depends both on the rescaled mass 
and the rescaled chemical potential, 



Veff 



2 _ 



mY^effV 



"eff 



(3.11) 



(3.12) 



We will also need the derivative of the condensate to generate the expectation value ((Tr[[/o + J/d])^)f/„. Whether 
or not to consider the ©(e^) correction is trivially achieved by keeping or dropping the subscript "e//" in the 
couplings in Eqs. p.lOp and p.l3p . For example, the partition function to leading order is given by Eq. p.3p 
with leading order couplings E and F. 

We define a second functional, which is the derivative of the partition function with respect to the rescaled 
chemical potential, to generate the following average: 



yAVeff,aeff) 



d 



ln[Z,] = (Tr[rC/trt/o])^ 



Uo 



We will also need a third generating functional 

MVeff,aeff) EE (Tr[(mo^mo)^])l^„ , 
to be given below in Eq. (j3.22|) in order to complete the calculation. 

3.1 Unflavoured correlation functions 

We begin with the scalar two-point function, inserting Eqs. (j3.7p - (j3.9p into Eq. (j2.26p to obtain 



(3.13) 



(3.14) 



'eff 



2S^ 



S2 



NfyA^,a) - (Tr[r])2 + — (A',(r;,a) - Nf) 



(3.15) 



Here 



. denotes the expectation value with respect to the full field U{x) at fixed topology ly. In the term 



proportional to A(x), the subscripts eff can be dropped since A(x) is already of 0{e^). The form of the first 
three terms is identical to the result for /i = in [1]. The only difference is that here, the resolvent £,^(77, a) also 
depends on the rescaled chemical potential. The remaining term proportional to is an explicit correction to 
the result of [3|, at non- vanishing chemical potential. 

For the pseudoscalar two-point function, we obtain in the same way 



Veff 



A{x) 



2E2 



dr] S 



Veff 



N, 



E2 r? E 



2 2 

+ ^{NfyAn,a) - (Tr[r])2 + '^{XAri^a) - Nj)) 



+ Nf 



+ Oie^) . (3.16) 



As before, all terms but the last are of the same form as for fi = 0, but now depend on the rescaled chemical 
potential a and thus F, also at leading order. The term in the last line proportional to is again an explicit 
correction term. 

We now explicitly give the three functions appearing in the results for the two-point correlations above in 
the particular case of Nf = 2 flavours. The results for more flavours follow easily using Appendix IbI 
The Nf = 2 flavour partition function reads [ini HO] 



Nf 



ZAVeff)^ / dAAe^"'//(4^'-2)/,(Ar;e//)\ 



(3.17) 
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where Ii, denotes the modified /-Bessel function. 

From the above equations, the condensate of Eq. p.lOp easily follows, and for two flavours we arrive at 



N 



f 



^^{r]eff,aeff) 



1 



ZuiVeff) Jo 



dXX 6 2 



*"'^^^^^'"''/.(Arye//)(— /.(A77e//) + A/,+i(A77,//)). (3.18) 



The second quantity 3^,y(??,a) equally follows from its definition p.l3p by differentiation. Using Eq. p.l7p we 
obtain for two flavours 



A^/ = 2 : y^iveff,aeff) 



ZuiVeff) Jo 



dXX^{AX' - 2)e^"=//(4^'-2)/.(A,?e//)' , (3.19) 



given here to C(e^), which we will need for the flavoured correlators below. 

The third function is obtained by differentiation of the following generalised partition function computed in 
Appendix |B] 



2„ 



J rr ^ .^rrv^ 3^f^f//M'Tr[r[/+r[/o]+iySe//™Tr[;7o+c/t] + c^Tr 
(IhUo det[(7g \e 



(3.20) 



It is remarkable that this group integral can be calculated, even when adding any power TA{TUlTUof\ to the 

partition function of Eq. p.3p , as was pointed out in [21] . We generalise the result given in [H] in Appendix |B] 
to arbitrary masses M. ^ nilNf, with Tr[r] ^ 0. For two flavours this generalised partition function reads 



Jo 



and for more flavours the result is given by Eq. (jB.Sp . Wc then obtain 

d 



(Tr {TU},rUof 



(3.21) 



(3.22) 



making only the ordinary partition function of Eq. (|3.17p appear in the denominator. Wc again display the 
two-flavour result at equal mass. 



^/ = 2 : XuiVeff,aeff) 



1 



ZuiVeff) Jo 



dXX (16(A4 - A^) + 2) e^''"}i^^^"-^h,{XT^,fff 



(3.23) 



3.2 Flavoured correlation functions 

We now give the sums over flavoured two-point correlations, which are easier to compute on the lattice, since 
they do not involve disconnected contributions. This is just a matter of inserting Eqs. (|3.7p - (|3.9p into the 
results of Eqs. (|2.26p and (|2.27p at fixed topology. Here it is useful to observe that we have already computed 

the combinations (Tr {Uo ± f/g)^ )uo ~ \(Tr[[/o ± ^o])^) ™ the unflavoured cases. We obtain for the 

scalar two-point function, 



Y^{Sa{x)Samv 



U(x) 



y2 



dVeff 



S2 



Veff 



"in 



+Nf + -^[NfyAVelf,(^eff) - (Tr[r])2 + --^{X,{T^,jf,a,fj) - Nj) 



2 , "e//, 



2F2 



+^/ + ]^(^/^- - (Tr[r])2 + y (A',(77,«) - Nf)) 



f S2 ' NfT] 

„2 



[N] - 2) 



77- 
(3.24) 
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and for the pscudoscalar correlator, 



^(Pa(x)Pa(0))^(,) 



^eff 



[N] - 1) i:^{r]^ff,aeff) 



Veff 
2 



a 



^eff 



A(x) 



2F2 



2 

{Njy,{ij,ff,a,jf) - (Tr[r])2 + ^ (A", a^// ) - Nf)) 

2 

(iV/:V.(r7,«) - (Tr[r])2 + y (A',(77,«) - Nj)) 



7V/772 



(3.25) 



where we now have explicit /z-dependent corrections at leading order, which is not the case for the other 
correlation functions. 



4 Conclusions 

We have computed scalar and pseudoscalar two-point correlation functions in the epsilon regime of Chiral 
Perturbation Theory, up to and including corrections of 0{t^). The main new feature of our results is the 
inclusion of an isospin chemical potential of real or imaginary type. This leads to the appearance of the low 
energy constant F already at leading 0{e^). We find corrections to the meson correlation functions obtained at 
/X = 0, which contain both explicit terms proportional to = fi^F'^V and a*, as well as implicit corrections 
which arise through those to the partition function and its derivatives. All of our results are obtained for 
A'^y-flavour QCD. Possible extensions of our work include the calculation of correlations functions in quenched 
or partially quenched QCD, which would require computing supersymmetric extensions of the group integrals. 
Axial or vector current two-point functions are also feasible. 

Our results provide alternate means of extracting the low energy constants S and F from lattice calculations. 
In particular, the additional parameter ^ and the fact that F appears at leading order should help improve the 
precision in the determination of this low energy constant. Moreover, we have obtained a new one-loop result 
for Feff, which includes one of the 0{e^) corrections to the partition function. 
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A Derivation of zero-mode group integral identities for /i 7^ 

In this appendix we will derive a set of U{Nf) group identities among expectation values of averages over various 
traces. Those relations are needed to arrive at Eqs. p.7p - p.9p in the expressions for the two-point functions. 
In particular we have to specify here 

C = /iP = /i diag(l7Vi,-l7V2) , (A.l) 

with A^i 7^ N2 in general. Furthermore let us stress that our identities hold for degenerate masses only, 

VJ:M = rj^Nf (A.2) 
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where 77 is the rescaled mass. In order to simphfy notation we will use 



O) = 



di/[/det[t/''] cxp 



ia2Tr[rC/tr[/] + i?7Tr[C/ + C/t] 



\- [ dnU O dctlU"] exp ]-a^T:T[TU'^TU] + ir/Tr[C/ + U'^]] 
-■f Ju{Nf) [2 2 



(A.3) 
(A.4) 



dropping most indices from the main text. In particular here U = Uq denotes the constant U{Nf) matrix. The 
full C(e^) improved expectation values Eqs. p.3p and (|3.4p are trivially obtained by shifting a — > a^ff and 

V Veff- 

Following [4] we introduce the explicit representation of the left differentiation with respect to group elements 
Uki of U{Nf) 

Va = i{taU)kl^ ■ (A.5) 



dUi 



kl 



Throughout this appendix, the ta denote the generators of the algebra of u{Nf). They satisfy the U{Nf) 
completeness relation 

{ta)ij{ta)kl = 2^il^jk 



in the normalisation Tr[iaib] = ■^Sab- For example, we obtain 

VaU = +itaU , VaC/t = -iUHa , 

Vadet[f7] = iTT[ta]Aci[U] . 
Due to the left invariance of the Haar measure, integrals over total derivatives with respect to Va vanish 



= 



dnU VaTr 



U{Nf) 



taG{U) AetlU"] exp 



(A.6) 
(A.7) 

(A.8) 



for any function G{U). By choosing a suitable set of functions, we will generate a closed set of equations for 
averages that can be solved for in terms of known generating functions. 

The simplest choice in Eq. (jA.Sp is G{U) — 1 and a = 0, without summing over a, which leads to the 
following identity 

= vNf + '^iTTp-U^) . 

This equation is invariant under complex conjugation when simultaneously changing v 
in Eq. (|A.4p ). Next, in choosing G{U) = U — U'' and keeping a = 0, we obtain: 



(A.9) 

-1/ (see the measure 



= |((Tr[C/ - U^f^ + uNf{TT[U - C/t]) + (Tr[C/ + C/t]) 
Next, we sum over a in Eq. (jA.Sp is G{U) — U, leading to: 



{Nf + ^)(Tr[C/]) + ll{Tr U' )-NA+-( (T 



U^TU'T 



{Tt[U]) 



Either by complex conjugation and changing —v, or simply choosing GiU) = —W we obtain 



= (7V^._^)(Tr[[/t]) + |((Tr[c/t2])-7V/)+^((Tr 



(Tr[[/t]) 



(A.IO) 



(A.ll) 



(A.12) 



Taking the sum of Eqs. (jA.ll|l and (|A.12p we obtain an equation for (Tr {U — L/^)^ 



The quantity (Tr 



and its conjugate which are new for a ^ cannot be easily derived from 

a known generating functional. Therefore we need an additional set of equations compared to [1], which is 
generated by choosing G{U) = UTU'T 



= 



(Tr 



) - (Tr[C/t])) + (Nf + z.)(Tr[[/r[/tr]) - ((Tr[r])2) + ^ ((Tr [([/rC/tr)^] ) - Nf) , 



(A.13) 



11 



and its complex conjugate 



{TT[U])j+{Nf-iy){TT[UTU^T])-^{TiiT]fj + — ((Tr [(C/FC/trf 



N, 



(A.14) 



We can now eliminate (Tr 



and its conjugate by using the sum of Eqs. (|A.13p and (|A.14|) to obtain 



(Tr[([/-C/t)2j) = _^(Tr[[/ + C/t])-— (Tr[[/-C/^]) 
L J ?7 ?; 



2a2 r 



2iV/(Tr[[/rC/tr]) -2^(Tr[r])2^ + a2(^(Tr[(C/r[/tr)2])-Ar^)] . (A.15) 



All objects on the right hand side can now be generated. Differentiating the partition function with respect to 
the mass we have 



{Tr[U + U^]) = 2|-ln[2.] = 2iV/^i^ 
using the definition Eq. (|3.10|) . Eq. (jA.Op provides the difference 



Differentiating the partition function with respect to we obtain (see the definition Eq. (|3.13p ) 



(A.16) 



(A.17) 



(A.18) 



The trace of the square is obtained from the generalised partition function Zg,,„ derived in Appendix IB] below. 

d 



(Tr {UTWT 



tF\2 



9. 



gen_ 



(A.19) 



see the definition Eq. ^J^. Inserting aU these into Eq. (jA.lSp we arrive at Eq. ([g^ . where (Tr ([/ + U^Y ) 
is obtained trivially by adding +4iV/. 

The missing squares of traces in Eqs. (|2.23p and (|2.25p follow by differentiating twice with respect to the 
mass, 



(Tr[t/ + C/t])^) . 4i-gz. . 4A.,A^^+47V?^^^^'")' 



(A.20) 



leading to Eq. p.7p . Finally, a last equation is needed that contains the square of the trace of the difference, 
which appears in Eq. (jA.lOp above. Using Eq. (jA.9p and the generator of Eq. (|A.16p . we obtain the Eq. ([3] 
in Sec. [21 after replacing the couplings a and rj by their effective counterparts, aeff and ?7e//- 



B Calculation of the generalised partition function 

In this appendix we compute the following generalisation of the partition function of Eq. (|3.20p 
-2^seri({?7}) = / dHf7odet[;7o]''exp 

JU{Nf) 

where = diag(77^i=i, • • • , j/atj , ?7/2=i, • ■ • i^ywa) contains the rescaled masses which may now be different. As 
in the previous appendix, the inclusion of effective couplings is trivial. Without loss of generality we choose 
N2 > Ni in Nf = A^i + A''2. The volume V and higher order coupling constants are all absorbed into the 
coefficients Up, where the sum can have a finite or an infinite number of terms. We only require that the 
integrals converge. 

The result for this generalised partition function was given in [21j for degenerate masses and A'^i = N2. Our 
generalisation follows [20] closely. Because we differentiate the generalised partition function with respect to 



-Tr[X((7o + U^)] + apTr \iUoTU^T)P 



(B.l) 
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the couplings ap in order to generate expectation values, we have to keep track of all constants that depend on 
the Up. The unitary matrix Uq can be parametrised as follows [TOl 



Uo 



A = 



vi 

V2 

A 



ui 

U2 



A 



- A2 



V 



In, - y 








vl 
vl 



-1n2-Ni I 



(B.2) 



Here, we denote by the matrix A = diag(Ai, . . . , A^Vi) containing the real numbers A^ G [0, 1] for k ~ 1, . . . , A^i. 
The unitary submatrices are G [/(A^i), U2 e U{N2) and vi G lJ(Ni) = U{N2)/{U{1)^^ x [/(TVs - A^i))- 

The matrix A is Hermitean and we observe that 



Tr 



[UoTulvy 



= Tr 



so that all unitary degrees of freedom drop out. From 

/ 

{KTf = 



2A2 



2MI1n, - A2 




(ArAr)f 



-2A^lAr, -A2 

2A2 - 1^, 




(B.3) 



(B.4) 



I7V2-W1 / 



we obtain 



Tr 



Tr 



2A2-1 



2A^ljVi -A2 

The 2Ni eigenvalues can be written in diagonal matrix form as 




Tr[ljV2-ArJ 



2A^ 



±i2X\ 1 



-Ni 



A2 . 



(B.5) 



(B.6) 



We thus arrive at 



Tr 



Ni 



Y.^T2p{X,) 



(B.7) 



The real polynomials T2p{Xi) of degree 2p that we obtain are the Chebyshev polynomials of the first kincfl]. 

Coming back to the integral in Eq. (jB.l|) . we can now go to an eigenvalue basis using the parametrisation of 
Eq. (|B.2p . Because of the decoupling of the unitary degrees of freedom in Eq. (|B.3p , the calculation is identical 
to the one presented in Sec. 3, to which we refer the reader for details. In particular integrating out all 
unitary degrees of freedom cancels the Jacobian from the parametrisation eq. (jB.2[) . Collecting all of these 
results, we obtain the following expression for Nf flavours: 



.(M) 



N 



AiVi({?7/i})A^,({77/2}) 



Z. 



gen (?7/l=l, '7/2=1 ) 



-■gen (Wi: '7/2=1) 
^!^('//2=l) 
'7/2=1 [Vf2=l) 



Zgen{-nfl=l,riN2) 
Zgen iVNi , 'nN2 ) 

Il'{VN2) 

v^r''^~'ii'''-'''~'\vN2) 



(B.8 



^ In 1211 the polynomial was given in the form cos (2pcos-l(A)). 
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where A^r^ — YijyiiVj ~Vi) is the Vandermondc determinant withm each flavour. The generahsed Nf 

flavour partition function, Zgen{vijV2), that is used inside the determinant above is given by 



dX X cxp 



^ap(2r2j,(A)-(7V2-A^i)) 



I^{X7]i)I„{Xr]2) 



(B.9) 



With I^\'rifi) we denote the fc-th derivative of the /-Bcssel function. The /-Bessel function itself, in fact, is a 
one-flavour partition function, Zi, ~ up to a trivial dependent constant. The constant N in Eq. (jB.Sp 

is an irrelevant normalisation factor that does not contribute in expectation values. Eq. (jB.Sp is the main new 
result of this appendix. It is the solution of the generalised partition function given by the group integral (|B.ip . 
It extends the previous result of [H] to nondegenerate masses and Tr[r] ^ 0. 
As an example, for p = 1, we get 

2r2(A) = 4A2 - 2 , (B.IO) 

which gives back the standard partition function of Eq. (|3.17|1 . for ai = \oi\^^ (and mass ?7e//)- It is explicitly 
given in Eq. (|3.17p . up to a constant prefactor. For p = 2, the second polynomial is 

2Ti{X) = 16(A*-A2) + 2, (B.ll) 

leading to the generalised partition function with a2 = oj that is needed in Eq. (I3.2ip . 
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